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This paper is on further development of discrete complex analysis introduced by R. Isaacs, R. Duffin, 
and C. Mercat. We consider a graph lying in the complex plane and having quadrilateral faces. A 
function on the vertices is called discrete analytic, if for each face the difference quotients along the 
two diagonals are equal. 

We prove that the Dirichlet boundary value problem for the real part of a discrete analytic function 
has a unique solution. In the case when each face has orthogonal diagonals we prove that this solution 
converges to a harmonic function in the scaling limit (under certain regularity assumptions). This 
solves a problem of S. Smirnov from 2010. This was proved earlier by R. Courant-K. Friedrichs- 
H. Lewy for square lattices, by D. Chelkak-S. Smirnov and implicitly by P.G. Ciarlet-P.-A. Raviart 
for rhombic lattices. 

In particular, our result implies uniform convergence of finite element method on Delauney tri- 
angulations. This solves a problem of A. Bobenko from 2011. The methodology is based on energy 
estimates inspired by alternating-current networks theory. 
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1. Introduction 

Various constructions of complex analysis on planar graphs were introduced by Isaacs, Duffin, 
Mercat [18, 15, 13, 22, 23], Dynnikov-Novikov [14], Bobenko-Mercat-Suris [2], Bobenko-Pinkall- 
Springborn [4]. Recently this subject is developed extensively due to applications to statistical 
physics [28], numerical analysis [17, 3], computer graphics [1, 29], combinatorial geometry [24]; see 
[21, 28] for recent surveys. 

This paper concerns linear complex analysis on quadrilateral lattices [2]. A quadrilateral lattice is 
a graph QcC with rectilinear edges such that each bounded face is a quadrilateral (not necessarily 
convex). Depending of the shape of faces, one speaks about square, rhombic, kite, or orthogonal 
lattices (the latter are quadrilateral lattices such that the diagonals of each face are orthogonal); see 
Figure 1. Different types of lattices are required for different applications; see Section 5. 
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Figure 1: Examples of lattices Q: (a) square; (b) rhombic; (c) orthogonal; (d) quadrilateral. 
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A complex- valued function / on the vertices of Q is called discrete analytic [23], if the difference 
quotients along the two diagonals of each face are equal, i. e., 

f(zi) ~ f(z 3 ) _ f(z 2 ) - f(z 4 ) 

Z\ — Z 3 Z2 — Z A 

for each quadrilateral face Z1Z2Z3Z4, see Figure Id. The motivation for this definition is that both 
sides of equation (1) approximate the derivative of an analytic function / inside this face. The real 
part of a discrete analytic function is called a discrete harmonic function. 

Discrete complex analysis is analogous to the classical complex analysis in many aspects [21]. One 
of the most natural and at the same time challenging problems is to prove convergence of discrete 
theory to the continuous one when the lattice becomes finer and finer [28]. A natural formalization 
of such convergence is uniform convergence of the solution of the Dirichlet boundary value problem 
for a discrete harmonic function to a harmonic function in the scaling limit. 

1.1. Previous work 

Convergence in this sense was proved by R. Courant-K. Friedrichs-H. Lewy [11, §4] for square 
lattices, by D. Chelkak-S. Smirnov [8, Proposition 3.3] and implicitly by P.G. Ciarlet-P.-A. Raviart 
[10, Theorem 2] for rhombic lattices. In fact convergence for rhombic lattices is equivalent to conver- 
gence of the classical finite element method [13]. The latter subject is well-developed; see a survey [7] 
and a textbook [6]. Nonrhombic lattices cannot be accessed by known methods. Weaker convergence 
results not involving boundary value problems were obtained in [22, Theorem 3], [17, Theorem 2]. 

1.2. Contributions 

We prove that the Dirichlet boundary value problem for a discrete harmonic function on a quadri- 
lateral lattice has a unique solution. Our main result is that in the case of orthogonal lattices this 
solution converges to a harmonic function in the scaling limit; see Convergence Theorem 1.2 below. 
This solves a problem of S. Smirnov [28, Question 1]. In concert we get a simpler proof for the 
particular cases known before. 

In particular, our main result implies uniform convergence of finite element method on Delauney 
triangulations; see Corollary 5.1 below. This solves a problem of A. Bobenko (private communication; 
see also [29, Table in Section 2]). 

1.3. Statements 

Let us give precise statements of main results. 

The boundary dQ of the graph Q is the boundary of its outer face. Hereafter assume for simplicity 
that dQ is a closed curve without self-intersections. Denote by Q° the set of vertices of the graph Q. 

Let g: C — > R be a smooth function. The Dirichlet (boundary value) problem on Q is to find 
a discrete harmonic function uq i9 : Q° — > R such that uq^(z) = g(z) for each vertex z G dQ. The 
function uq >9 : Q° — > R is called a solution of the Dirichlet problem. 

Uniqueness Theorem 1.1. The Dirichlet boundary value problem on any finite quadrilateral lattice 
has a unique solution. 

This theorem is nontrivial because discrete harmonic functions do not satisfy the maximum 
principle in general; see Example 3.6 below. 

Let Q C C be a domain. The Dirichlet (boundary value) problem on Q is to find a continuous 
function uq^: CIQ — > R harmonic in Q and such that un }9 (z) = g(z) for each point z e dQ. The 
harmonic function un t9 : Q — > R is called a solution of the Dirichlet problem. 

A sequence of lattices {Q n } approximates the domain Q, if for n — > 00: 

• the maximal distance from a point of dQ n to the curve dQ tends to zero; 

• the maximal edge length of Q n tends to zero. 
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A sequence of lattices {Q n } is nondegenerate uniform, if there is a constant Const (not depending 
on n) such that for each member of the sequence: 

(D) the ratio of the diagonals of each face is less than Const and the angle between them is greater 
than 1/ Const; 

(U) the number of vertices in an arbitrary disk of radius equal to the maximal edge length is less 
than Const. 

A sequence of functions u n : — > C converges to a function u: Q — > C uniformly on each 
compact subset, if for each compact set K C Q we have max zgA - n QO \u n (z) — u(z) \ — > as n — > oo. 

Convergence Theorem 1.2. Let Q C C be a domain bounded by a smooth closed curve dfl without 
self-intersections. Let g: C — > R be a smooth function. Let {Q n } be a nondegenerate uniform 
sequence of finite orthogonal lattices approximating the domain Q. Then the solution UQ n , g : Q° n — > R 
of the Dirichlet problem on Q n converges to the solution u^, g '■ ^ — > R of the Dirichlet problem in Q 
uniformly on each compact subset ofQ. 

1.4- Organization of the paper 

In Section 2 we introduce main ideas of the proofs and state key lemmas. In Sections 3 and 4 
we prove Uniqueness Theorem 1.1 and Convergence Theorem 1.2, respectively. In Section 5 we give 
applications of our results to numerical analysis, network theory, probability theory, and state some 
open problems. 

2. Main ideas 

2.1. Energy minimization 

Our approach is based on energy estimates inspired by alternating-current networks theory. 
Recall that the (Dirichlet) energy of a continuous piecewise-smooth function it: Q — > R is 

E n (u) := / \Vu\ 2 dA. 
Jn 

This is a convex functional on the space of continuous piecewise-smooth functions with fixed boundary 
values, and harmonic functions are characterized as minimizers of this functional. 

Let us define a discrete counterpart of the energy, which is the main concept of the paper. The 
gradient of a function u: Q° — > R at a face zi £2-23^4 of the quadrilateral lattice Q is the unique vector 
Vqu(ziZ2Z 3 Z4) G R 2 such that 

V Q u(z 1 z 2 z 3 z 4 ) ■ z^z% = u(zi) - u(z 3 ), 
Vqu(ziz 2 z 3 z a ) ■ Z2~z\ = u(z 2 ) - u(z A ). 

The energy of the function u : Q° — > R is the number 

E(u) := ^ l^Quiz^z^)] 2 ■ Area(z 1 z 2 z 3 z 4 ), (2) 

21Z22324CQ 

where the sum is over all the faces 2:1^2^3^4 of the lattice Q. 

We give a physical motivation for this definition in Section 5.2. A similar but nonequivalent 
definition was given in [23, Formula (12)]. Our energy has the same properties as its continuous 
counterpart: 

Convexity Principle 2.1. The energy E(u) is a strictly convex functional on the affine space 
MP ~ 9 ® of functions u: Q° — > R having fixed values at the boundary dQ. 

Variational Principle 2.2. A function u: Q° — > R has minimal energy E(u) among all the functions 
with the same boundary values if and only if it is discrete harmonic. 

These principles are proved in Section 3. Uniqueness Theorem 1.1 is their direct consequence. 
After the "right" discrete energy has been guessed, these results are proved by standard methods. 
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2.2. Energy estimates 

Let introduce more delicate energy estimates required for the proof of Convergence Theorem 1.2. 

Joining the opposite vertices in each quadrilateral face of the lattice Q, we get two connected 
graphs B and W associated to the lattice; see Figure 2. (The vertices are joined by a straight line 
segment, if the segment lies inside the face, and by the 2-segment broken line through the midpoint 
of the opposite diagonal, otherwise.) The eccentricity of a lattice Q is the infimum of the numbers 
Const such that the lattice satisfies conditions (D) and (U) from Section 1.3. Throughout the paper 
we use the following notation: 

• Q is an arbitrary domain bounded by a smooth curve without self-intersections; 

• g : C — > R is an arbitrary smooth function; 

• B and W are the two graphs associated to the lattice Q. 

• e is the eccentricity of the lattice Q; 

• h is twice the maximal edge length of the lattice Q. 




Figure 2: The graphs B and W associated to a quadrilateral lattice Q. 



Energy Convergence Lemma 2.3. Let {Q n } be a nondegenerate uniform sequence of quadrilateral 
lattices approximating the domain f2. Then E(g |qo ) — > E^{g) as n — > oo. 

Equicontinuity Lemma 2.4. Let Q be an orthogonal lattice, K be a compact set inside dQ, and 
u: Q° R be a discrete harmonic function. Denote by r := min 2e ^ Dist(z, dQ). Then there is a 
constant Const x, r ,e depending only on K, r, e (but not on Q, u, z, w) such that 



\u(z) — u(w) 



< Const^ r e • E(u) 1/2 ■ ln~ 1/2 (1 + , T , ) 

\ \z — w\J 



(3) 



for any z, w G K D B°. 

The laplacian Aqu: Q° — > R of a function u: Q° — > R is defined by the formula 

dE{u) 



du{z) 



for each z G Q°. 



Laplacian Approximation Lemma 2.5. Let Q be a quadrilateral lattice and R be a square of side 
length r > h inside dQ. Then there is a constant Const e depending only on e (but not on Q, g, R, 
r, h) such that 



[A Q (g\ Q o)](z)- f AgdA 

-DnDIl JR 



z£RnB° 



< Const e ( hrmax\D g(z) \ + r max \D g(z) 

z£R zeR 
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For a subset K G C and a function u: B° — > K denote -£^(w) := XLei^nB () w^- 2 )- 

Friedrichs Inequality Lemma 2.6. Le£ Q be a quadrilateral lattice, K be a compact set inside dQ, 
and u: B° — » K. be an arbitrary function. Denote by r := max zg ^ Dist(^, dQ). Assume that r > h 
and r > max 2e gQ Dist(z, dQ). Then there is a constant Constn je depending only on Q and e (but not 
on Q, K, r, h, u) such that 

h 2 L 2 n _ K (u) < Const Qje (hrL 2 dQ (u) +r 2 E(u)) . 

These four results are proved in Section 4. Lemmas 2.3, 2.5, 2.6 are proved using suitable 
modifications of the approaches of [12], [8], [11], respectively. Equicontinuity Lemma 2.4 is essentially 
novel. Estimates analogous to (3) were known for square lattices [11, equation (12) in §4.2], for 
rhombic ones [8, Corollary 2.9, Proposition 2.7 and Appendix A], and for orthogonal ones [25, 
Corollary 3.4]. The methods known for square and rhombic lattices do not generalize to more 
general ones (because in general there are no discrete exponentials [19] and no higher derivatives of 
discrete analytic functions). Surprisingly, our proof is simpler than the ones in [8, 11] even for the 
particular types of lattices studied there. 

Sketch of the proof of Convergence Theorem 1.2 modulo the above lemmas. (See the details in Sec- 
tion 4.7.) Take a compact set K dQ. Restrict each function uq u ^ to the set B® f]K. By Variational 
Principle 2.2 and Energy Convergence Lemma 2.3 we have E(uQ ntg ) < E(g \qo) — > E$i(g) < oo for 
n — t- oo. Thus the sequence E(uQ n , g ) is bounded. Then by Equicontinuity Lemma 2.4 it follows 
that the sequence of functions uq u ^ is equicontinuous. By the Arzela-Ascoli theorem a subsequence 
UQ k ^ g of the sequence mq„ i9 converges uniformly to some function u: Q — > R. Using the Weyl lemma 
and Laplacian Approximation Lemma 2.5 we show that the function u is harmonic. Using Friedrichs 
Inequality Lemma 2.6 we estimate J n _ K {u — g) 2 dA and thus verify the boundary condition. We 
conclude that u equals the unique solution of the Dirichlet problem. Thus the initial sequence 
u Q n ,g ( n °t just a subsequence) converges to uniformly on each compact subset of Q. □ 

3. Uniqueness 

In this section we prove Uniqueness Theorem 1.1 and the results stated in Section 2.1. We also 
prove two results (Maximum Principle 3.5 and the Green Identity 3.7) required for the next section. 

3.1. Convexity Principle 

Proof of Convexity Principle 2.1. Consider the linear space MP° of functions u: Q° — > M,. Let the 
coordinates of a function u : Q° — > R be the values of the function. Clearly, the gradient Vw linearly 
depends on u and thus the energy E{u) is a quadratic form in u. So it suffices to prove the convexity 
of E{u) in the case when the affine space MP C MP passes through the origin, that is, all the 
fixed boundary values equal zero. 

Clearly, E(u) > for each u e R Q 9Q . It remains to prove that E(u) = only if u — 0. Assume 
that E(u) = 0. Then Vqu(ziZ2ZsZ4) = for each face ^1^2-23-24 of the lattice Q. This means that 
for each face Z\ Z2Z3Z4 we have u(z\) = u(z^) and u(z2) = u^z^). Any face can be joined with the 
boundary dQ by a sequence of faces such that two neighboring ones share a common edge. Thus for 
each vertex z e Q° there is a boundary vertex w G dQ such that u(z) = u{w). Since the boundary 
values equal zero it follows that u = 0. □ 

3.2. Variational Principle 

Denote by * : M 2 — > M. 2 the counterclockwise rotation through 7r/2 around the origin. Two 
functions u, v : Q° — > M. are conjugate, if Vqv = * Vqu. 

Claim 3.1. Two functions u,v : Q° — > M. are conjugate if and only if the function u + iv : Q° — > C 
is discrete analytic. 



5 



Proof. Identify the gradient VqU G M 2 with a complex number Vqu G C. The function u + iv: Q° — > 
C is discrete analytic if and only if 

u(zi) + iv(zi) - u(z 3 ) - iv(z 3 ) _ u(z 2 ) + iv(z 2 ) - u{z±) - iv(z 4 ) 

Z\ — Z 3 Z 2 — ^4 

for each face z\z 2 z 3 z± of the lattice Q. Substitute the expression 



u{zi) - u(z 3 ) = VqU ■ (z x - Z 3 )/2 + VqU • {zx ~ Z 3 )/2 

and analogous ones for v(zi) — v (z 3 ), u(z 2 ) — u^z^), u(z 2 ) — u(z 4 ) into the above equation. We get 



V z\ — z 3 z 2 — z^j 

Since the second factor in the left-hand side is nonzero, the equation is equivalent to Vq1> = iVqu. □ 
Claim 3.2. A function u: Q° — > M. has a conjugate if and only if for each z G Q° — dQ we have 

^2 *V Q u{z x z 2 z 3 z±) ■ z^z% = 0, 

21222:324 : 21=2 

where the sum is over all the faces z\z 2 z 3 z 4 of the lattice Q such that z\ = z. 

Proof. Let us prove the "only if" part. Assume that v. Q° — > K is conjugate to u. Then 

21222324:21=2 21222324:21=2 21222324:21=2 

because the diagonals of type £2-24 form a closed cycle around a nonboundary vertex z. 

Let us prove the "if" part. Denote by V(z 2 z^) := *Vqu ■ z±z\. Assume that for each bounded 
face W\W 2 . . .w m , w m+ i := Wi, of the graph W we have J2i<k<m^( w k w k+i) — 0. Then a function 
v : W° — > R is well-defined by the formula v{w m ) := Y^i<k<m V{ w k w k+i)i where Wiw 2 ...w m is a 
path in the graph W joining w m with a fixed vertex ui\. Define a function v. B° — > R analogously. 
Consider the combined function v. Q° — > 1R. Then for each face Z\z 2 z 3 z± of the lattice Q we have 
v{zi) — v{z 2 ) = V{z 2 z±) = *Vqu-z^z\ and analogously v{z 3 ) —v{zi) = *Vqu- z 3 z\. Thus Vqv = *Vqu, 
and v is conjugate to m. □ 

Claim 3.3. For each vertex z G Q° we have 

[A Q u\(z) = ^2 *V Q u{ziz 2 z 3 z 4 ) ■ zjzt. 

21222324 : 21=2 

Proof. Let v : Q° — > R be equal to 1 at the vertex z and at all the other vertices. Differentiating 
the energy E{u) and applying the identity (c • d)(*a ■ b) = (*c • b)(a ■ d) — (*c ■ a){b ■ d) we get 

[A Q u)(z)= ^2 {V Q u ■ V Q v)(-2Axea > (z 1 z 2 z 3 z i )) 
21222324CQ 

21222324CQ 

= ^2 ((*^Q u ' z i z s){VqV ■ ~ (*V q m • z^ztj^QV ■ zTzt)) 
21222324CQ 

= ^2 {{*Vqu-z^z 3 ){v{z 2 ) -v{z a )) - {*V Q u ■ z^zt)(v(zi) - v(z 3 ))) 

21222324CQ 

= 2J *VqM • Z^z\. 

21222324 : 21=2 

□ 

Proof of Variational Principle 2.2. By Claims 3.1-3.3 a function is discrete harmonic if and only if 
its laplacian vanishes at nonboundary vertices. By Convexity Principle 2.1 the latter is equivalent 
to having minimal energy among the functions with the same boundary values. □ 
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3.3. Uniqueness Theorem 

Proof of Uniqueness Theorem 1.1. By Convexity Principle 2.1 the energy E: MP — » K has a 
unique global minimum u G M9 ~ 9 Q. By Variational Principle 2.2 the function w is the solution of 
the Dirichlet problem and it is unique. □ 

Remark 3.4. Define a discrete Riemann surface to be a cell decomposition Q of a surface with 
quadrilateral faces together with an identification of each face with a quadrilateral Z\z 2 z^Zi C C by 
an orientation preserving homeomorphism. (No agreement of such identifications for different faces 
is assumed.) The results of Sections 2.1 and 3 remain true for an arbitrary simply-connected discrete 
Riemann surface, not necessarily a quadrilateral lattice in the complex plane. 

3.4- Maximum Principle 

Let us discuss the case of orthogonal lattices in more detail. For a face z\z 2 z^z^ with the vertices 
listed clockwise denote by c{z\z?) := In the case of an orthogonal lattice we have c{z\z$) > 0, 

and the energy and the laplacian take the usual form 

E{u) = ^2 M 2 ^) -«0i)) 2 + c(z 1 z 3 y 1 (11(34) -u{z 2 )f) /2; 

21 22 23 24 

[Aqu](zi) = c(zi* 3 ) (u(z 3 ) - w(zi)) • 

23 : 2123CB 

Thus the value of a discrete harmonic function u at a nonboundary vertex of 5 equals to the weighted 
mean of the values at the neighbors. This immediately implies the following known result. 

Maximum Principle 3.5. Let Q be an orthogonal lattice and let u: Q° — )■ K be a discrete harmonic 
function. Then 

maxn(z) = max u(w) and maxu(z) = max u(w). 

z£Q° w£Q°ndQ 2g_B° w£B°ndQ 

For an orthogonal lattice Uniqueness Theorem 1.1 is an immediate corollary of Maximum Prin- 
ciple 3.5 and the finite-dimensional Fredholm alternative. However, the principle does not hold for 
nonorthogonal lattices in general. 

Example 3.6. (S. Tikhomirov) Let M > 1 and let Q be the lattice formed by the quadrilateral with 
the vertices 0, tan(7r/8), M(tan(-7r/8) + and the 3 other quadrilaterals obtained by symmetries 
with respect to the origin and the coordinate axes; see Figure 3. Define a discrete analytic function 
/ : Q° ->■ C by the following table: 








±z ±tanf ±M(tanf + z) ±M(tanf-t) 


/(*) 




10 Mi 



Set u := Re /. Then maxgo uj max^g u = M/2 and thus can be arbitrarily large for large M. If B 
is the graph formed by 4 diagonals from the origin then max^o u > whereas max^o n gQ u = 0. 

Mi\^^\^^*0 

O.^^^^^^A.Mi 

Figure 3: A discrete analytic function on a nonorthogonal quadrilateral lattice not satisfying the maximum principle. 
The values of the function are shown near the vertices. 
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3.5. The Green Identity 

Let us state one more result specific for orthogonal lattices, which is required for the sequel. 

Green Identity 3.7. Let Q be an orthogonal lattice andu,v: B° — > C be arbitrary functions. Then 



z€B° 

Proof. For an orthogonal lattice the energy splits as E(u) = Eb(u) + Ew{u), where Eb(u) and 
Ew{u) depend only on the values of the function u at the vertices of B and W, respectively. For an 

arbitrary homogeneous quadratic form Eb{u) we have J2zeB° ( u ( z ) d dv{z) ~ v ( z ) 9 du{z) ^) = ^' w hich 
is equivalent to the required identity. □ 

For nonorthogonal lattices the Green identity does not remain true unless one replaces summation 
over B° by summation over Q°. 

4. Convergence 

In this section we prove Convergence Theorem 1.2 and the results stated in Section 2.2. 

4-1. Geometric preliminaries 

Let us start with some basic estimates involving the lattice eccentricity. For a subgraph R C B 
denote by E R {u) the sum (2) over all the faces Z\Z 2 z 3 z 4 of Q containing an edge of R. 

Path Energy Claim 4.1. For any path wqW\ . . . w m C B we have 

(u(w m ) -U(W )) 2 

Proof. Estimating |Vqw| through its projection, using condition (D) from Section 1.3, the inequality 
sin(l/e) > l/2e, and the Schwarz inequality we get 

E W0 ... Wm (u) = \ VQu\ 2 Are&(z 1 z 2 z 3 Z4) 

^ (u(w k ) - u(w k -i)) 2 i I, \ • /( \ 
> ) : 7^ Zi2 3 \z 2 z 4 \smZ(z 1 z 3 ,Z2Z 4 ) 



k=l 



\ZlZ 3 \ 



>A (u(w k ) - u{w k _i)) 2 . 1 

> > ■ sin - 

^ e e 

k=l 

(u(w m ) - u(wq)) 2 

~ 2me 2 

where the first sum is over all faces z\Z2Z 3 z 4 of Q containing an edge of the path Wq . . . w m . □ 
Projection Claim 4.2. For any face ziz 2 z 3 z 4 of the lattice Q and any vector it we have 

\1T- 



~J\ <4e ' , , 31 +4e 



ziz\\ \lf ■ z 2 z\\ 



\ZlZ 3 \ \Z2Z 4 \ 



Proof. Since Z(ziZ 2 , z 3 z\) = ±Z( v , z 3 z\) ± Z(17, z 3 z 4 ) it follows that at least one of the angles in 
the right-hand side, say, the first one, does not belong to the interval (| — 7L, | + ^-)- Then 



1 


v ■ ZtZ 3 \ 


2e 


Z\Z 3 \ 





< 4e 



v ■ ZlZz\ 

\ziz 3 \ 

□ 



8 



Rectangle Capacity Claim 4.3. A rectangle r x h with the side r > h contains at most Aer/h 
vertices of the graph B. 

Proof. The rectangle r x h can be covered by 4[r//i] discs of radius h/2. Then by condition (U) from 
Section 1.3 the number of vertices in the rectangle is less than 4e[r/h}. □ 

Diameter Claim 4.4. The diameter of each bounded face of the graphs B and W is at most h. 

Proof. A bounded face of the graph B contains a vertex of the graph W. The vertex is joined by 
edges of the graph Q with all the vertices of the face (and by "half-edges" of the graph W with the 
break points of the 2-segment edges of B in nonconvex faces of Q). Since the edges of Q have length 
at most h/2 it follows that the diameter of the face of B is at most h. □ 

4-2. Convergence of energy 

For the proof of Energy Convergence Lemma 2.3 we need the following claim. 

Gradient Approximation Claim 4.5. We have \Vg — Vq(# \qo)\ < 8eh max zg conv(aQ) \D 2 g(z)\. 

Proof. Consider a face z^z^z^ of the lattice Q. By the Rolle theorem there is a point z G 2:12:3 
(possibly outside z\ -22-23^4 but inside the convex hull Conv(<9Q)) such that (Vg(z) — [Vq^] {z\ Z2Z3Z4)) ■ 
zizt/\zizt\ = 0. Thus (Vg - V Q g) • z^/\z^\ < /imax 2eConv ( zlZ2Z324) \D 2 g{z)\ in the face z 1 z 2 z 3 z 4 . 
The same inequality holds with Z\Z 3 replaced by -22-24. By Projection Claim 4.2 the claim follows. □ 

Proof of Energy Convergence Lemma 2.3. Denote by Q the domain enclosed by the curve dQ. Since 
Q n approximates Q it follows that some neighborhood Q' of Q contains all the lattices Q n and 
Area(f2 — Q n ), Are&(Q n — Q) — > as n — > 00. Since the domain Q is bounded and the function 
g: C — > R is smooth it follows that Vg is bounded in Conv(fi'). Thus the integrals Eci(g), Efi n (g) 
exist and En(g) — E^ n (g) = E Q _Q n (g) — Eg n _ n (g) — > as n — > 00. By Gradient Approximation 
Claim 4.5 we get E^^g) — E(g |qo ) — > as n — > 00, and the lemma follows. □ 

Remark 4.6. For a harmonic function v. Q — > R and a sequence of functions v n : Q — > R conditions 
v n — v on dfl and En(v n ) —> Eq(v) do not necessarily imply that v n — > v pointwise. For instance, 
take a continuous function v n : D 2 — > R such that v n (z) = for \z\ = 1, v n (z) = 1 for \z\ < l/n, 
and v n (z) is harmonic in the ring l/n < \z\ < 1. Then E^(v n ) = 1/lnn — > as n — > 00 but 
u n (0) = 1^0 = v(0). 

4-3. Equicontinuity 

Proof of Equicontinuity Lemma 2.4 for square lattices. First assume that B is a square lattice and 
the segment joining z and w is contained in the graph B D K; see Figure 4. Assume without loss of 
generality that u(z) > u(w). 

For now denote by h the step of the square lattice B. Let R m be the boundary of the rectangle 
2mh x (2mh + \z — w\) centered at the point (z + w)/2 with the side 2mh orthogonal to the segment 
zw; see Figure 4. 

















h 
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Figure 4: The rectangles R m and the points z m ,w m from the proof of Lemma 2.4. 
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Take m < r/2h. Then the distance from R m to the segment zw C K is less than r/y/2. So R m 
lies inside dQ and thus R m C -B. Then by Maximum Principle 3.5 there are points z m ,w m G R m 
such that u(z m ) > u(z) and u{w m ) < u{w). The points z m and u> m are joined in the graph R m by a 
path of length at most Am + |z — w\/h. Since the eccentricity of a square lattice is 1, by Path Energy 
Claim 4.1 we get the following estimate for the energy on the subgraph R m 

. . > K^j -u(w m )\ 2 > |u(g) -m(w)| 2 _ 
fiml J ~ 8m + 2|2 - iw|//i ~ 8mA + 2\z - w\ ' 

Summing these inequalities for m from to [r/2h] and estimating the sum via an integral we get 

[r/2h] 1 / r \ 

E(u) > E RM > -\u(z) - u(w)\ 2 \n 1 + ^— . 

This is equivalent to the required inequality. □ 

For nonsquare lattices the proof is essentially the same but requires more technical details. 

Proof of Equicontinuity Lemma 2.4 for convex compact sets. Let us prove the lemma in the case 
when the segment zw C K. 

Consider an auxiliary square lattice with edges parallel and orthogonal to zw and with step equal 
to the maximal edge length h of the lattice Q. Define the rectangles R m literally as above. 

For each m = 1,2, ... , [r/2h] the boundaries of the rectangles R m and R m -i are separated by 
a simple closed path R m C B. Indeed, otherwise i? m and R m -\ would be joined by a path in 
the complement to the graph B and thus the graph B would have a face of diameter > h, which 
contradicts to Diameter Claim 4.4. 

Number of edges in the path R m is not greater than the number of vertices of the graph B lying 
in the strip between R m and R m +i- Thus by Rectangle Capacity Claim 4.3 the path R m contains 
less than 4e(m + \z — w\/h) edges. Now the same energy estimates as for the square lattice (with 
R m replaced by the path R m ) prove the required inequality. □ 

Proof of Lemma 2.4 in the general case. For each point of K take a disc of radius r/2 around the 
point. The interiors of these discs form an open covering of K. Since K is compact it follows that 
the covering contains a finite subcovering. Let Di, . . . , Dn be the closures of the discs of the finite 
subcovering. Denote K' := U^ =1 Dk- Then K C K' and K' is inside dQ. 

Assume without loss of generality that dK' does not have "cusps", i.e., the discs Di, . . . , D N do 
not touch each other. Assume also that K' is connected. Denote by Dist(z, w) the distance between 
points z and w in K', i. e., infimum of lengths of all broken lines joining z with w and lying in K' . 

Let us prove the asymptotic form Dist(^,io) < Const/^ r |z — w\ for some number Const^,. de- 
pending on K and r. Indeed, denote by e the minimal distance between two disjoint discs among 
Di, . . . , Dn- li\z — w\ < e then the points z and w are contained in a pair of intersecting discs Di and 
Dj (or possibly in the same disc). Since Di and Dj do not touch each other the required asymptotic 
form follows. For \z — w\ > e denote <p(z,w) := Dist(z, w)/\z — w\. Since is a continuous function 
on a compact set { (z, w) G K' x K' : \z — w\ > e }, the required asymptotic form follows again. 

Now take two points z, w G K' . Join them by a broken line of length less than Dist(£, w) + \z — w\. 
Since K' is a union of N discs we may assume that the broken line has at most iV edges, and each 
edge is contained in one of the discs. Since Dist(z,w) = Const^ r (|z — w\) it follows that we can 
subdivide the broken line to get a broken line z = z±, z 2 , ■ ■ ■ zm+i = w with M < Const Kr + 1 edges 
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of length < \z — w\ each. Then 

M 

\u(z) - u(w)\ < ^ \ u i z k) ~ u(z k+1 )\ 

k=l 

M 

< Const e ■ E{ufl 2 ln- 1/2 (l + r\z k - z k+1 \- 1 ) 

k=l 

< Const^, r , e • E{u) l/2 lrT 1/2 (l + r\z - w]- 1 ). 

Here the first inequality is obvious. The second asymptotic form follows from the particular case of the 
lemma proved before because each edge z k Zk+i is contained in one of the discs D\, . . . , D N . The third 
asymptotic form holds because the function ln _1 ^ 2 (l + r|z| _1 ) is increasing and M < Const jc, r +l- d 

Remark 4.7. Our proof of Equicontinuity Lemma 2.4 cannot be generalized to nonorthogonal lattices 
because it essentially uses Maximum Principle 3.5 not holding for them. 

4-4- The Friedrichs inequality 

For the proof of Friedrichs Inequality Lemma 2.6 we need some notation and auxiliary claims. 
Let abed be a rectangle with vertices a,b,c,d G C listed clockwise such that ab is outside dQ and 
\b — c| > h; see Figure 5. Denote r := \b — c\. For a point z G abed denote by R z the rectangle 
(r + 2h) x 2h centered at (V + w')/2 with the side 2h parallel to ab, where z' G cd and w' G ab are 
the points such that z G z'w' and z'w' \\ be. 




Figure 5: The rectangle abed and the neighborhood R 



Claim 4.8. A vertex z G -B fl abed can be joined with dQ by a path in the graph R Z C\ B. 

Proof. Assume that there is no path as required. Then the point z is separated from dQ by a path 
P in the complement R z — B. Since ab is outside dQ it follows that z is separated from the point 
w' by the path P as well. The path P is not closed because the graph B is connected. Thus the 
endpoints of the path P belong to dR z . The path P intersects the segment zw' because P separates 
z from w'. Thus the face of the graph B containing P has diameter greater than h, a contradiction 
to Diameter Claim 4.4. □ 

Claim 4.9. We have w G R z if and only if z G R w . 

Proof. Straightforward. □ 
Claim 4.10. h 2 L 2 abcd (u) < 8erhL 2 dQ (u) + 64e 4 r 2 £(V). 
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Proof. Take a vertex z G abed D -B . By Claim 4.8 it follows that the vertex z is joined with a vertex 
w G dQ by a path in the graph R z (IB. Take such path with minimal number m of vertices. Then m 
is not greater than the total number of vertices of the graph B in the rectangle R z . Hence m < 4er/h 
by Rectangle Capacity Claim 4.3. Thus by Path Energy Claim 4.1 we get 

u{z) 2 < 2u{w) 2 + 2(u(z) - u(w)) 2 < 2L 2 9QnRz (u) + 16e 3 rE BnRz (u) /h. 

Sum these inequalities over all z G abed fl B°. A vertex w G dQ can contribute to L 2 d Q nR ^(u) only if 
w G R z . Thus by Claim 4.9 and Rectangle Capacity Claim 4.3 it contributes at most 4er/h times. 
Similarly, an edge z"w" of the graph B can contribute to Ebhr z (u) only if z" G R z . Thus by Claim 4.9 
and Rectangle Capacity Claim 4.3 it contributes at most 4er/h times. Thus our summation leads to 
the required inequality. □ 

Proof of Friedrichs Inequality Lemma 2. 6. Since dQ is a smooth curve it follows that its r-neighborhood 
can be covered by finitely many rectangles abed such that ab is outside the r-neighborhood of Q and 
\b — c| = 3r. Moreover, the number of rectangles is bounded by a number depending only on the do- 
main Q (but not on K). These rectangles cover the strip Q — K as well. Since r > max zg gQ Dist(z, dQ) 
and r > h it follows that for each rectangle the side ab is outside dQ and \b — c\ > h. Summing the 
inequalities of Claim 4.10 for each of the rectangles, we get the required inequality. □ 

Remark 4.11. Our proof of Lemma 2.6 does not generalize to domains with nonsmooth boundaries, 
e.g., to a domain bounded by a cardioid. 

4-5. Approximation of laplacian 

First we prove Laplacian Approximation Lemma 2.5 in several particular cases and then combine 
them together. Throughout this subsection z denotes the coordinate in the complex plane C, so 
that, e.g., Re z denotes the function z^Kez. 

Proof of Lemma 2.5 for g(z) = Rez. Clearly, the function f(z) = z is discrete analytic. Thus g(z) = 
Kez is discrete harmonic. Then by Variational Principle 2.2 it follows that [Aq Hez](w) = for each 
vertex w G Q° — dQ. On the other hand, A Re z = as well, and the lemma follows. □ 

Proof of Lemma 2.5 for g(z) = Imz. This is proved by the previous argument with iz instead of 
z. □ 

Proof of Lemma 2.5 for g(z) = \z\ 2 . For a face zxz^z^z^ of the lattice Q denote by z' the intersection 
point of the middle perpendiculars to the diagonals z\z^ and z^z^. Clearly, then Vq\z—z'\ 2 [zxz^z^ z±) = 
0. Using the two previous cases of the lemma and Claims 3.1, 3.3 we get 

[A Q \z\ 2 }(w) = [A Q \z - w\ 2 ](w) + [A Q (M 2 - 2Rewz)}(w) 
= (*V Q \z-w\ 

21222324 : 21 =U) 



Z4Z2 



(2*V Q Re((z'-w)z) + *V Q \z 

21222324 : z±=w 

2 Vq lm{{z'-w)z) ■ z^zt 

2324 : 21=10 

2Im ({z' -w){ Zi - z 2 ) S j 

2324 : z\=w 

= 4Are&(wz 2 z' Z4). 



z'\ 2 1 • 2:42:2 



212:22:324 : 21=10 



21222324 : zi=w 



212:22324 : z\=w 



Here and in the next paragraph the area of a closed broken line is understood in oriented sense. 
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Denote by Rq the closed broken line formed by the edges z 2 z' and z'z^ for all the faces ziz^z^ 
of Q such that z\ G R and z 3 (jL R. By Projection Claim 4.2 it follows that z' Z\ < eh, thus Rq is 
contained in the e/i-neighborhood of the curve dR. Summing up the above expressions over all the 
vertices w G B° D R we get 



[A Q \z\ 2 ](w)- f A\z\ 2 dA 



weB°nR. 



4Area(i? Q ) - 4Area(i?) < Voehr. 

□ 



For the next case of the lemma we need the following two claims. 

Claim 4.12. For an arbitrary function u: W° — > C we have ^2 zlZ3 ,bQji(u(z 2 ) ~ u ( z 4}) — 0, where the 
sum is over all the faces Z1Z2Z3Z4 such that z\ G B° fl R, z 3 G" R. 

Proof. Take an arbitrary vertex z 2 G W° — dQ. It is contained in a face w 1 . . . w m of the graph 
B. Assume that the vertices of the face are listed in clockwise order. Let us move along the path 
w\ . . . w m . Each time we cross the boundary dR the value u(z 2 ) appears in the considered sum with 
certain sign. The sign is positive, if we move from the exterior of R to the interior at the moment, 
and negative — otherwise. Since both dR and W\ . . . w m are closed curves it follows that the vertex 
z 2 appears in the considered sum equal number of times with positive sign and with negative sign, 
and the claim follows. □ 

Claim 4.13. The number of faces ziz 2 z 3 Z4 such that z\ G B° (1 R, z 3 (jL R is at most 96er/h. 

Proof. By Diameter Claim 4.4 it follows that the vertices of such faces are contained in the h- 
neighborhood of dR. By Rectangle Capacity Claim 4.3 the number of these vertices is at most 
32er//i. Then by the Euler formula for planar graphs it follows that the number of faces is at most 
9Qer/h. □ 

Proof of Lemma 2.5 for g(z) = Re z 2 . Analogously to the previous case of the lemma we get 

[A Q Re z 2 }{w) = 21m((z'-w) (z±-z 2 )). 

Summing up this expressions over all w G B°f]R, canceling repeating terms and applying Claim 4.12 

we get 

[& Q Rez 2 }(w)= 2lm(z'(z 4 -z 2 )) = Yl ^((^' - z 2 - z A )(z 4 - z 2 )) . 

w&RnB zxz 3 \\\dR ZiZ 3 (hdR 

Here the first sum is over all the vertices z\ G B° PI R and the other sums are over the faces z\Z 2 z 3 z A 
of Q such that z\ G R and z 3 G" R. By Projection Claim 4.2 it follows that \z' — z 2 \, \z' — z±\ < eh. 
Thus by Claim 4.13 we get 



J2 [A Q Rez 2 ]( Zl )- j ARez 2 dA 

- DO J R 



ieR.nB a 



< Y \2z'-z 2 -z 4 \-\zi-z 2 \ < 96er/h-2eh-h < 192e 2 hr. 

ziz^fridR 



□ 

Proof of Lemma 2.5 for g(z) = Imz 2 . This is analogous to the previous case. □ 
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Proof of Lemma 2.5 in the case when D k g = at the center of R for k = 0, 1, 2. Using the estimate 
< 4rmax ze ij \D 3 g(z)\ we get 



Ag dA 



R 



< 4r 3 max \D 3 g(z)\. 

z£R 



Now applying Claim 3.3, canceling repeating terms, applying Claim 4.13, Gradient Approximation 
Claim 4.5, and the estimate |Vg(z)| < r 2 m&x zeR \D 3 g(z)\ we get 



z!eRnB° 



E E *^ Q9 ' ^ 



21222324 : 21=2 



E *^Q9 ■ ^2 

2l23(tl9i? 

< 96er//i ■ (r 2 + Aehr) max |-D 3 _g r (^) | ■ ft 

zeR 

< 480e 3 r 3 max |L> 3 #(z)|. 

2gi? 



□ 



Proof of Lemma 2.5 in the general case. Assume without loss of generality that the center of the 
square R is the origin. By the Taylor formula 

g(z) = a + Oi Re 2; + a 2 Imz + a 3 Re^ 2 + a 4 |2;| 2 + a 5 Imz 2 + 

for some a 1; . . . , a 5 G R and a function such that D k g(0) = for k = 0, 1, 2. For each of the 
summands the lemma has already been proved. Thus the general case follows. □ 

4-6. Uniform limit 

The following lemma is the last result we need before the proof of Convergence Theorem 1.2. We 
write x n ~ y n and x n >z y n , if lim(x n — y n ) = and lim(x n — y n ) > 0, respectively. 

Lemma 4.14. Suppose that under the assumptions of Convergence Theorem 1.2 the sequence UQ n:9 |#o 
converges to a continuous function u: Q — > R uniformly on each compact subset ofQ. Then u = Mn,g- 

To prove the lemma, first let us establish harmonicity of the limit function. 

Claim 4.15. Let {Q n } be a nondegenerate uniform sequence of orthogonal lattices approximating 
a domain Q. Let u n : Q° n — )■ R be a sequence of discrete harmonic functions. Suppose that the 



restrictions u n | #o converge to a continuous function u : Vt — > 
Q; then the function u is harmonic. 



uniformly on each compact subset of 



Proof of Claim 4-15. Take an arbitrary smooth function v : Vt — > R vanishing outside a compact 
subset K C Q. By the Weyl lemma it suffices to prove that J n uAvdA = 0. Assume without loss of 
generality that dK is smooth and v = also in a neighborhood of dK. 

Let us estimate the difference between J Q uAvdA and its discrete counterpart. For each n take 
an auxiliary infinite square lattice with edge length r := \/~h. For a face R of the n-th auxiliary 
lattice denote u n (R) := ma.x z€RnK u(z). Then u n =4 u on the compact set K because u is continuous 
on this set. Applying the convergence u n ,u n =4 u, the boundness of Av,Aq n v, and then Laplacian 
Approximation Lemma 2.5 we get 



I uAvdA- ^[u n A Qn v]{ 



R : RDK^tl) 



I AvdA- E t A Q^K 

Jr zeRnBO 



< 



Area(K)/r 2 ■ max \u\ ■ Const e • (rh + r 3 ) max \D 2 v , D 3 v\ < Constg^^^ • \/h — > as n — > 00. 

K K 
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It remains to estimate the discrete counterpart of J n uAv dA. Take n large enough so that K is 
inside dQ n . Applying Green Identity 3.7 and the assumptions that u n is discrete harmonic and v 
vanishes outside K, we get 

Thus j n uAv dA = 0, which proves the claim. □ 

Proof of Lemma 4-14- By Claim 4.15 the function u is harmonic. It remains to verify the boundary 
condition. 

Denote by S(r) C Q is the set of points with distance to the boundary dQ less than r. First 
let us prove the estimate h 2 L 2 s ^(uQ nt9 — g) < Const ej Q ig r 2 for r > h and some number Const ei n iS 
not depending on r. Apply Friedrichs Inequality Lemma 2.6 for the function mq„ i9 — g. The energy 
E(uQ nt9 —g) < 2E(uQ n: g)+2E(g qo ) is bounded by Variational Principle 2.2 and Energy Convergence 
Lemma 2.3. We have I} d Q n {uQ n ^—g) = because mq„ iS = g at dQ n , and the required estimate follows. 

Now by Diameter Claim 4.4 it follows that for each p < r 

/ (u - gfdA ■< h 2 L 2 s{r) _ s{p) (u Qn:9 - g) < h 2 L 2 s{r) (u Qn:9 - g) < Const eA9 r 2 . 

JS(r)-S(p) 

Approaching p — > we get £ fs( r )( u ~ 9) 2 dA — > as r — > 0. By [11, §4.1] this condition implies the 
boundary condition. □ 

4-7. Convergence Theorem 

Proof of Theorem 1.2. Take an arbitrary subsequence Q nk of the given sequence of lattices Q n . For 
brevity denote Qk '■= Qn k , B k := B nk . Take a sequence of compact sets K 1 C K 2 C • - • C ft such 
that Vt = Uf =l Kj. 

Let us estimate \uq^ 9 \. Since the sequence Qk approximates the domain Q it follows that there 
is a disk Q' containing all lattices Qk- By Maximum Principle 3.5 we have 

m ^\uQ k , g ( z )\ = m n a x \uQ h ,g( w )\ < \g(w)\ < oo 

because UQ kt9 = g at dQk and (7: C —¥ R is continuous. So the sequence MQ fcj5 is uniformly bounded. 

Let us estimate the right-hand side of inequality (3) from Equicontinuity Lemma 2.4 for K :— K\ 
and u := UQ k>g . Since the sequence Qk approximates the domain Q it follows that there is a number 
k\ such that for each k > k\ the set K\ is contained inside dQk and Dist(iTi, dQk) is bounded 
from zero. By Variational Principle 2.2 and Energy Convergence Lemma 2.3 we have E(uQ k g ) < 

E(g qo ) — > En(g) < 00 for k — )■ 00. Thus the sequence E{uQ k ^ g ) is bounded. 

Then by Equicontinuity Lemma 2.4 it follows that UQ kig 

positive function 5(e) not depending on k and such that for each z,w G K fl with \z — w\ < 5(e) 
and k > k\ we have \uQ kt9 (z) — UQ k}9 (w)\ < e. By the Arzela-Ascoli theorem it follows that there is a 
function u\ : K\ — > C and a subsequence {Ik} of the sequence 1,2,... such that k\ is as before and 
uq u9 converges to u\ uniformly in K\. 

Now proceed to the next compact set K 2 . Analogously, there is a function w 2 : K 2 — > C and 
a subsequence {m^} of the sequence {4} such that mi — l±, m 2 — l 2 , and MQ mi9 converges to u 2 
uniformly on K 2 . Clearly, u\ = u 2 on K\. Thus the extension can be continued, and eventually we 
get a function u: Q — > C and a subsequence {pk} of the sequence 1,2,... such that UQ pt9 converges 
to u uniformly on each compact subset of Q. 

By Lemma 4.14 the function u is the solution u^^ of the Dirichlet problem in Q. Since the 
solution uq j9 is unique it follows that the initial sequence UQ nt9 : — > M converges to u^ >9 uniformly 
on each compact subset of Q. Analogously, UQ ni9 : W° — > R converges to uq !9 uniformly on compact 
sets. This completes the proof of main results. □ 



3 o is equicontinuous, i. e., there is a 
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5. Applications and open problems 



5.1. Application to numerical analysis 

Convergence Theorem 1.2 provides a new approximation algorithm for numerical solution of the 
Dirichlet boundary value problem. It also gives a new convergence result for the classical finite 
element method, which we are going to state now. 

Let B be a triangulation of a polygon B. The finite element method approximates the solution 
of the Dirichlet problem on Q by a continuous function ub, 9 '- B — > R which is linear on each face 
of B, equal to the given function g on the boundary B° D dB, and has minimal energy E B (u) = 
J^\Vu\ 2 dA (among such functions). Equivalently [12, §4], the restriction ub, 9 - B° — > R can be 
defined as the unique function equal to g on B° R dB and such that for each z\ G B° — dB we have 
J2 ZS c ( z iZ3)(uB,g(zi) —UB, g {z^)) = 0, where the sum is over all neighbors z 3 of z\, and c{z\Z 3 ) is given 
by the formula 

c(z 1 z 3 ) := (cot a + cot (5) /2, (4) 

where a and (3 are the angles opposite to the edge ziz 3 in the two triangles of B sharing the edge; 
see Figure 6. The function m^ s : B° — > R is called the solution of the Dirichlet problem on B. 







B / 













Figure 6: A kite lattice Q associated to a Delauney triangulation B. 

Usually one proves convergence of the finite element method under certain assumptions on in- 
dividual triangles [7]. For instance, it was proved in [9, Theorem 3.3.7] that us n ,g convergences 
uniformly to u% g1 if there is a constant Const such that 

(A) the minimal angle of each triangle is greater than 1 /Const; 

(R) the ratio of any two edges of each triangulation is less than Const. 

According to [7] no uniform convergence results without assumptions (A) and (R) were available. 

Following [5] we suggest a new approach measuring "triangulation quality" via configuration of 
neighboring triangles rather than the shape of individual ones. A triangulation B is called Delauney, 
if a + (3 < 7r for each pair of adjacent triangular faces (and thus c(ziz 3 ) > above). A Delauney trian- 
gulation exists for any prescribed set of vertices B° [5]. A sequence of triangulations is nondegenerate 
uniform, if there is a constant Const such that for each member of the sequence 

(D) for each edge the sum of opposite angles in the (two or one) triangles containing the edge is 
less than 7r — 1/Const (in particular, the triangulation is Delauney); 

(U) the number of vertices in an arbitrary disk of radius equal to the maximal edge length is less 
than Const. 

Assumption (U) is weaker than (R); neither (D) nor (A) is weaker than the other one. We prove 
convergence of the finite element method for triangulations satisfying (D) and (U): 
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Corollary 5.1. Let Q C C be a domain bounded by a smooth closed curve dQ without self- 
intersections and let u: C — > R be a smooth function. Let {B n } be a nondegenerate uniform sequence 
of triangulations approximating the domain Q. Then the solution ub„, 9 '- B® — > R of the Dirichlet 
problem on B n converges to the solution : fi — > R of the Dirichlet problem on Q uniformly on 
each compact subset ofQ. 

Proof. In each triangle of B n , draw 3 segments joining the circumcenter with the vertices. Erase 
hanging edges from the obtained graph. We get an orthogonal quadrilateral lattice Q n . Formula (4) 
gives the same values c(z\Z 3 ) as in Section 3.4; see Figure 6. Thus UB n)9 '- B® — > R is the restriction 
of a discrete harmonic function Q® t — > R. By Uniqueness Theorem 1.1 we have ub^ 9 {z) = UQ n:9 (z) 
for each z G B^. Since {B n } is a nondegenerate uniform sequence of triangulations approximating 
the domain Q, the same is true for {Q n }. By Convergence Theorem 1.2 the corollary follows. □ 

Vice verse, the finite element method can be applied to establish convergence of discrete harmonic 
functions on quadrilateral lattices. Using standard finite element described above one can approach 
only rhombic lattices (and also kite ones at the cost of establishing convergence only at the vertices 
of the graph B but not W). 

The following nonconforming finite element might be useful in the case of general quadrilateral 
lattices. Given a function u: Q° — > R define its interpolation Iqu: Z\Z 2 z 3 Z\ — > R to be the linear 
function on a face z\z%z 3 z± of Q such that [Iqu](zi) = u(zi), [Iqu](z 3 ) = u(z 3 ), and [Iqu](z2) — 
[Iqv\{z^) = u(z 2 ) —u(z4). Combining such linear functions together we get a (discontinuous) function 
Iqu: Q — )• R on the union Q of all the quadrilateral faces of Q. Clearly, then E{u) = E^(Iqu). 
Remark 5.2. In the case when all the bounded faces of Q are convex one can prove that E(u) = 
Eq(Ibu, Iwu), where Ibu,Iwu: Q — > R are piecewise linear extensions of the function u: Q° — > R 
to the faces of the graphs B U Q and W U Q, respectively. 

Problem 5.3. Give an effective approximation algorithm for finding the solution of the Dirichlet 
problem on a (nonorthogonal) quadrilateral lattice. 

Problem 5.4. Estimate the rate of convergence in Theorem 1.2. 

5.2. Physical interpretation 

Classical physical interpretation of complex analysis on orthogonal lattices uses direct-current 
networks [13] (for elementary introduction to networks see [26, 27, 24]). Let us give a new physical 
interpretation for arbitrary quadrilateral lattices involving alternating-current networks. This gives 
some insight and also interesting in itself. 

Define admittance of an edge C B by the formula 

c(z!Z 3 ) := i— — — , (5) 
z\ — Z3 

where Z1Z2Z3Z4 is the face containing Z1Z3 with the vertices listed clockwise. Clearly, this number has 
positive real part (and in case of an orthogonal lattice it is simply a positive number). 

A graph B with edge admittances having positive real parts can be considered as an alternating- 
current network; see [24, §2.4] and [16]. 

Given a discrete analytic function /: Q° — > C, define the voltage drop V(ziZz) and the current 
I(ziz 3 ) on an oriented edge z^zj, C B by the formula 

V(ziz 3 ) := f(zi) - f(z 3 ), I(z 1 z 3 ) := if(z 2 ) - if(z 4 ), 

where z±z 2 z 3 z 4 is the face of the lattice Q with the vertices listed clockwise. Boundary voltage 
drops are the differences f\zi) — f(z 3 ) for all pairs of consecutive boundary vertices zi,z 3 e B° D 
dQ. Boundary currents (or incoming currents) are the values if{zi) — if{z2) for all the pairs of 
consecutive boundary vertices z 2 ,z i e W° fl dQ. The voltage drop at a moment t is the number 
Re {V(ziz 3 ) exp(it)); the current at the moment t is defined analogously. 
A reformulation of Uniqueness Theorem 1.1 is the following result. 
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Corollary 5.5. Boundary voltage drops at the initial moment and boundary currents after one 
quarter of the period uniquely determine all the voltage drops and currents in an alternating- current 
network at all the moments of time. 

Problem 5.6. What is the physical meaning of the boundary condition, which is the Dirichlet one 
at the initial moment and the Neuman one after one quarter of the period? 

The physical meaning of Convergence Theorem 1.2 is that the voltage in a distributed direct- 
current network can be approximated by voltages in lumped direct-current networks. 

This physical interpretation gives also one more motivation for the definition of energy from 
Section 2.1. The energy of the network (dissipated per period) is 



where the sum is over all the edges ziz^ C B. Expressing the energy through u = Re/ and the 
admittances c(z\Zs), we arrive at the definition from in Section 2.1. 

5.3. Probabilistic interpretation 

Let Q be an orthogonal lattice and let B be one of the graphs obtained by joining the opposite 
points in each quadrilateral face of Q. A random walk on the vertices of B is defined as follows. At 
each moment of time the walker moves from his current position to one of the neighboring vertices 
with the probability proportional to the weights of the corresponding edges given by formula (5). 

The results of the present paper allow to generalize many estimates from [8] to nonrhombic 
lattices. Let us sketch two particular problems; see [8, 20] for accurate definitions. The probability 
that a random walk starting at a vertex z G B° first exits the domain Q through an arc A C dfl 
is called the discrete harmonic measure u)a.,si,q{z). A trajectory of a loop-erased random walk is 
obtained from a trajectory of a random walk by deleting loops in chronological order. 

Problem 5.7. Prove that for a nondegenerate uniform sequence of orthogonal lattices {Q n } approx- 
imating the domain Q (so that Q is inside each curve dQ n ) the discrete harmonic measure 0JA,si,Q n { z ) 
converges uniformly to its continuous counterpart. 

Problem 5.8. Prove that the trajectories of loop-erased random walks on orthogonal lattice converge 
to SLE 2 curves in the scaling limit. 

5.4- Generalizations 

Problem 5.9. Generalize Theorem 1.2 to: 

(1) nonorthogonal quadrilateral lattices; 

(2) nonuniform sequences, i.e., not satisfying condition (U) from Section 1.3 (for adaptive meshes); 

(3) singular boundary values (for convergence of discrete harmonic measure, the Green function, the 
Cauchy and the Poisson kernels, abelian integrals); 

(4) domains with rough boundaries (for probabilistic applications); 

(5) other types of boundary conditions; 

(6) other Riemannian surfaces; 

(7) other elliptic PDE. 

Problem 5.10. Prove that under the assumptions of Convergence Theorem 1.2 the gradient VQiiQ n)9 
converges to Vuq )9 uniformly on each compact subset of Q. 

Problem 5.11. Construct a sequence of quadrilateral lattices approximating a planar domain such 
that the solutions of the Dirichlet problem on the lattices do not converge uniformly to the solution 
of the Dirichlet problem in the domain. 

Problem 5.12. For which nonorthogonal lattices Q Maximum Principle 3.5 remains true? 




(6) 
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